QED calculation of the nuclear magnetic shielding for hydrogen-like ions 
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We report an ab initio calculation of the shielding of the nuclear magnetic moment by the bound 
electron in hydrogen-like ions. This investigation takes into account several effects that have not been 
calculated before (electron self-energy, vacuum polarization, nuclear magnetization distribution), 
thus bringing the theory to the point where further progress is impeded by the uncertainty due to 
nuclear-structure effects. The QED corrections are calculated to all orders in the nuclear binding 
strength parameter and, independently, to the leading order in the expansion in this parameter. 
The results obtained lay the ground for the high-precision determination of nuclear magnetic dipole 
moments from measurements of the <?-factor of hydrogen-like ions. 
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I. INTRODUCTION 

Last years showed much progress in the experimental 
determination of the ^-factors of hydrogen-like ions [l|- 
Measurements, accurate up to a few parts in 10 11 
, were performed by studying a single ion confined in a 
Penning trap. These experiments provided the stringent 
tests of bound-state quantum electrodynamics (QED) 
theory and yielded the best determination of the elec- 
tron mass @ ■ 

In order to match the experimental accuracy, many so- 
phisticated calculations were performed during the past 
ears, in particular those of the one-loop self-energy 0- 
the one- loop vacuum-polarization Q, the nuclear re- 
coil [nj, and the two-loop QED effects p], GJ]. These 
calculations made it possible to determine the electron 
mass from the experimental values of the bound-electron 
g factor [l3j]. The theoretical accuracy of the bound- 
electron g factor in hydrogen-like ions is presently at the 
1CP 11 level [12J for light elements up to carbon but dete- 
riorates quickly for heavier elements because of the un- 
known higher-order two-loop QED effects scaling with 
the nuclear charge number Z as Z . 

The experimental investigations of the ^-factors of 
hydrogen-like ions were so far performed for ions with 
spinlcss nuclei. However, when applied to the ions with 
a non-zero nuclear spin, such investigations can be useful 
as a new method for the determination of the magnetic 
dipole moments of the nuclei. This method has impor- 
tant advantages over the more traditional approaches, 
such as nuclear magnetic resonance (NMR), atomic beam 
magnetic resonance, collinear laser spectroscopy, and op- 
tical pumping (OP). These advantages are that (i) the 
simplicity of the system under investigation (a hydrogen- 
like ion) allows for an ab initio theoretical description 
with a reliable estimation of uncalculated effects and (ii) 
the influence of the nuclear-structure effects (which are 
the main limiting factors for theory) is relatively weak. 
This is in contrast to the existing methods, in which the 



experimental data should be corrected for several physi- 
cal effects, which are difficult to calculate. Among such 
effects is the diamagnetic shielding of the external mag- 
netic field by the electrons in the atom. The NMR results 
should be also corrected for the paramagnetic chemical 
shift caused by the chemical environment [l4j and the 
OP data are sensitive to the hyperfine mixing of the en- 
ergy levels [To| . Significant (and generally unknown) un- 
certainties of calculations of these effects often lead to 
ambiguities in the published values of nuclear magnetic 
moments |l6[. 

The goal of the present investigation is to perform an 
ab initio calculation of the <?-factor of a hydrogen-like ion 
with a non-zero nuclear spin. It can be demonstrated 
that the nuclear-spin dependent part of the atomic g- 
factor can be parameterized in terms of the nuclear mag- 
netic shielding constant a, which describes the effective 
reduction of the coupling of the nuclear magnetic mo- 
ment ft to an external magnetic field B caused by the 
shell electron(s), 
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The relativistic theory of the g-factor of a hydrogen- 
like ion with a non-zero nuclear spin (and, thus, the the- 
ory of the nuclear magnetic shielding) was examined in 
detail in Ref. [I?} • In the present work, we go beyond 
the relativistic description of the nuclear magnetic shield- 
ing and calculate the dominant corrections to it, namely, 
the self-energy, the vacuum-polarization, and the nuclear 
magnetization distribution corrections. As a result, we 
bring the theory to the point where the uncertainty due 
to nuclear-structure effects impedes further progress. 

The main challenge of the present work is the calcu- 
lation of the self-energy correction. To the best of our 
knowledge, the only previous attempt to address it was 
made in Ref. [l8| . In that work, the self-energy contri- 
bution to the shielding constant was estimated by the 
leading logarithm of its Za expansion (where a is the 
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fine-structure constant). In our work, we calculate the 
self-energy correction rigorously to all orders in the bind- 
ing nuclear strength parameter Za and, independently, 
perform an analytical calculation to the leading order in 
the expansion in this parameter (including both the log- 
arithm and constant terms). First results of this work 
were reported in Ref. |l9j . 

The rest of the paper is organized as follows. In Sec.llTl 
we summarize the relativistic theory of the g-factor of an 
ion with a non-zero nuclear spin and the theory of the 
nuclear magnetic shielding. Our calculation of the self- 
energy and vacuum-polarization corrections to all orders 
in Za is described in Sec. IIIII In Sec. |IV] we report the 
calculation of the QED correction to the nuclear mag- 
netic shielding to the leading order in the Za expansion. 
Sec. [V] deals with the other effects, namely, the nuclear 
magnetization distribution, the nuclear recoil, and the 
quadrupole interaction. Numerical results and discussion 
are given in Sec. IVII The paper ends with conclusion in 
Sec. EIIl 

The relativistic units (to = h = c = 1) and the charge 
units a = e 2 /(47r) are used throughout this paper. 



II. LEADING-ORDER MAGNETIC SHIELDING 

We consider a hydrogen-like ion with a non-zero spin 
nucleus placed in a weak homogenous magnetic field B 
directed along the z axis. Assuming that the energy 
shift due to the interaction with the B field (the Zee- 
man shift) is much smaller than the hyperfine-structure 
splitting (hfs), the energy shift can be expressed in terms 
of the g factor of the atomic system gp, 



AE = gp no B Mp , 



(2) 



where B = \B\, (1q — |e|/(2m) is the Bohr magneton, e 
and to are the elementary charge and the electron mass, 
respectively, and Mp is the z projection of the total an- 
gular momentum of the system F. To the leading order, 
the energy shift is given by 



AE = (FM F \ 



-—(fx a) 
2 v ; 
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where \FMf) = \jIFMp) is the wave function of the ion, 
j and / are the angular momentum quantum numbers of 
the electron and nucleus, respectively; /2 is the operator 
of the magnetic moment of the nucleus. The matrix ele- 
ment ([3]) can easily be evaluated, yielding the well-known 
leading-order relativistic result [21| . 
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where g^ is the Dirac bound-electron g factor [2Cj , gi = 
hKhnF) is the nuclear g factor, fx = (II\fl\II) is the 
nuclear magnetic moment, fi^ = \e\/(2m p ) is the nuclear 



magneton, m p is the proton mass, and (j ■ F) and (I ■ F) 
are the angular-momentum recoupling coefficients, 

(j-F) = [F(F + l)-I(I + l)+j(j + l)}/2, (5) 

(I-F) = [F(F + l)+I(I+l)-j(j + l)]/2. (6) 

Generalizing the leading-order result to include the 
higher-order correction, we write the atomic g factor gp 
as 



(j-F) to (I-F) 
^-^> (F + 1) -— 3/(1"^) F(F + 1) 



(7) 



In the above equation, the bound-electron g factor gj — 

gj +a/ (2ir) + . . . incorporates all corrections that do not 
depend on the nuclear spin, whereas the nuclear shield- 
ing constant a parameterizes the nuclear-spin dependent 
corrections. The bound-electron g factor gj has been ex- 
tensively studied during the last years, both theoretically 
0-[f3] and experimentally [H-Q- The goal of the present 
work is the ab initio theoretical description of the nuclear 
shielding parameter a. 

It can be seen from Eq. ([7]) that the nuclear-spin depen- 
dent contribution to the atomic g factor gp is suppressed 
by the electron-to-proton mass ratio and thus is by about 
three orders of magnitude smaller than the nuclear-spin 
independent part proportional to gj. It is, however, pos- 
sible to form a combination of the atomic g factors which 
is free of the nuclear-spin independent contributions. So, 
for ions with the nuclear spin / > 1/2, we introduce the 
sum of the atomic g factors g that is directly proportional 
to the nuclear magnetic moment, 



9 = 9f=i+i/2 + 9f=i-i/2 = -2 
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This combination of the g factors is particularly conve- 
nient for the determination of the nuclear magnetic dipole 
moments from experiment. Indeed, if both the gF=i+i/ 2 
and gp = j_i/ 2 <?-factors are measured and a is known from 
theory, Eq. ([8]) determines the nuclear magnetic moment 
fx. For the ions with a nuclear spin / = 1/2, Eq. (JSj) is not 
applicable and the nuclear magnetic moment should be 
determined from Eq. ([7]). 

Contributions to the nuclear magnetic shielding are de- 
scribed by Feynman diagrams with two external interac- 
tions, one with the external magnetic field (the Zeeman 
interaction), 



V zee (r)= l -^-B-(fxa) : 



(9) 



and another, with the magnetic dipole field of the nucleus 
(the hfs interaction) 



t r / \ e ^ rx a 

Vhfe(r) = -j-m 3- 
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(10) 



The leading-order contribution to the magnetic shielding 
comes from the following energy shift 

AE = 2 V ^^( a |^ zcc |,i)(n|y hfs | a ) , (11) 
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where the summation runs over the whole Dirac spec- 
trum with the reference state excluded. As follows from 
Eqs. ([2]) and 0, contributions to the shielding constant 
Sa are obtained from the corresponding energy shifts 8E 
by 
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For the electronic states with j a — 1/2, all nuclear 
quantum numbers in Eq. (fTTj) can be factorized out. The 
expression for the shielding constant is then obtained 
from that formula by using the following substitutions 

f r ,-, — » \ t r 7- r (fX Q?)o 

V zee V zcc = (r x a)o , V hfs -> V h{s = y ^ /u , 

|o) |oy 3 ) , |n) -> |ny 2 ) , 

2^a (prefactor) , (13) 

where the zero subscript refers to the zero spherical 
component of the vector. Here and in what follows, 
|m/ 2 ) = \n n ,fj, n = 1 /2) denotes the Dirac state with the 
rclativistic angular quantum number n n and the fixed 
momentum projection /i„ = 1/2. 

So, for j a = 1/2, the leading contribution to the mag- 
netic shielding is given by a simple expression 

er (0) =a V (ay 2 |\4cc|«i/ 2 ) {ny 2 |Vhfs|ai/ 2 ) . (14) 



Performing the angular integrations with help of 
Eq. (|A3|) . one gets for the ns reference state 



d(-2) 



(15) 



K n — — 1,2 n^a 

where are the radial integrals of the form 

poo 

Rab = / drr 2 + a [g a (r)f b (r) + f a (r)g b (r)] , (16) 
Jo 

g(r) and f(r) are the upper and lower radial components 
of the Dirac wave function, respectively, and the angular 
prefactors x Kn are given by x Kn —-i = —2/3, x K „=2 = 
-V2/3. 

For the point nuclear model, the sum over the Dirac 
spectrum in the above expression can be evaluated ana- 
lytically p3 - l24l | , see also more recent studies [13, HI] , 



r(°) 



4aZa ( 1 1 



9 V 3 6(1 + 7) 7 2 7 -l 

/ 1 97 
a Za - + —{Zaf + . 
U 108 v ' 



(17) 



where 7 = yl — (Za) 2 . For the extended nucleus, the 
calculation is easily performed numerically [17] . 



III. QED CORRECTION 

A. Self-energy: General formulas 

The self-energy correction in the presence of the ex- 
ternal magnetic field and the magnetic dipole field of 
the nucleus is graphically represented by six topologi- 
cally non-equivalent Feynman diagrams shown in Fig. [T] 
General expressions for these diagrams are conveniently 
obtained by using the two-time Green's function method 
[26l | . The resulting formulas are summarized below. 



1. Perturbed- orbital contribution 

The irreducible contributions of diagrams on 
Figs. [Ua)-(c) can be represented in terms of the 
one-loop self-energy operator E. This part will be 
termed as the perturbed-orbital contribution. The 
corresponding energy shift is 

AEpo = 2(E * (e a -H)' Vzcc (e a -HY Vhis) 

+2 fcpby Vhis T^h)> Vzcc) 

-2 (E R - 1 tt ^, V Ms }(V zcc ) 



(Sa ~ Hf 



-2(Z R )(V ZC 
-2 (V„ 



1 



(e„ - H) 



V S -R 7 — Vhfs) 



{e a -H)> n {e a -B)> 
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where we used the short-hand notations for the reduced 
Green's function, 



1 = \ n )( n \ 
(e„ - H)> L> 



and its derivative 



E 



\n)(n\ 



(e„ - H) 2 ' ^ (s a - e n ) 2 



(19) 



(20) 



The (unrenormalized) self-energy operator £(e) is de- 
fined by its matrix elements as follows 



(i|E(e)|fc) = ^ / du> 



(in\I (uj)\nk) 



£ — UJ — US r 



, (21) 



where I(lu) is the operator of the electron-electron in- 
teraction, I(lu) = e 2 a l _ l a u D tll '{uj), D^ v (uo) is the photon 
propagator, are the Dirac matrices, and u = 1 — iO, 
where iO is a small imaginary addition that defines the 



(a) (b) (c) (d) (e) (f) 

FIG. 1: Self-energy correction to the nuclear magnetic shielding. Double line represents the electron in the binding nuclear 
field. The wave line terminated by a triangle represents the dipole hyperfine interaction with the nucleus and the wave line 
terminated by a cross represents the interaction with the external magnetic field. 



positions of the poles of the electron propagators with re- 
spect to the integration contour. The summation over n 
runs over the complete Dirac spectrum. Renormalization 
of the one-loop self-energy operator is well described in 
the literature, see, e.g., Ref. [27]. In this work, the renor- 
malized part of the self-energy operator is denned as 

E R (s) = £(e) -/3<5ro- (e-a-p- V c -(3m)B {1) , 

(22) 

where 5m is the mass counterterm, j3 is the Dirac (3 ma- 
trix, is the one-loop renormalization constant, Vc 
is the binding Coulomb potential of the nucleus, and the 
renormalization is to be performed in momentum space 
with a covariant regularization of ultraviolet (UV) diver- 
gencies. Details on the renormalization procedure and 
explicit formulas for Sfl(e) can be found in Refs. [28l.[29|. 



where is the one-loop renormalization constant and 
the renormalization is to be performed in momentum 
space with a covariant regularization of UV divergencies, 
see Ref. [2!| for details. 

The irreducible contribution of the diagram shown 
in Fig. HJe), together with the corresponding derivative 
term, is referred to as the Zeeman-vertex contribution. It 
is given by 



AE V 



= 2(r 



i 
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{Sa - H) 



7 Vhfs 



+2 (£' 



R 



(e a - H) 



7 Vms) (Kec) , (26) 



where the 3-point Zeeman vertex is defined analogously 
to the hfs one. 



2. Single-vertex contributions 

The irreducible contribution of the diagram shown in 
Fig. [TJd), together with the corresponding derivative 
term, is referred to as the hfs-vertex contribution. It is 
given by 



.his 



2(r 



i 



hfs.-R 



■v„ 



+2 (E' fl — 1 — V zcc ) (VW , (23) 

where Fhf s ,_R = rhfs,_R(e a ) is the renormalized part of the 
3-point vertex representing the interaction with the hfs 
field and H' R = T,' R (e a ) is the derivative of the renor- 
malized self-energy operator over the energy argument, 
= ^/(^ £ )^fl( £ )l e = £ ■ The unrenormalized 3- 
point hfs vertex operator is defined by its matrix elements 
as 



i f°° 

(i\T h[s {e)\k) =— duo 

{in 2 \I{ui)\n 1 k)(n 1 \VM S \n 2 ) 



E 



(e - w - ue ni )(e — w - ue ri2 ) 



(24) 



The renormalized part of the operator is obtained as 

T h{s , R (e) = r hts (e) - V Ms , (25) 



3. Double-vertex contribution 

The contribution of the diagram shown on Fig. [TJf), 
together with the corresponding derivative terms, will be 
termed as the double-vertex contribution. It is defined by 



A£ d . vr - 2 (A(e„)) 

+ (X")(V zcc )(V his ) + (r> Ms )(V Zi 



1 " f 



aa'\I(ui)\a'a} 



(27) 



where A denotes the 4-point vertex representing the in- 
teraction with both the Zeeman and hfs interactions, 



(i|A(e)|fc) 



i 

2^ 



dco 



niTi2ii 3 v 



(m 3 |/(a;)|7iifc)(ni|V r ZC c|n2)(»-2|Mifs|»-3) 
-w — ue ni )(e — w — ue n2 )(e — u — ue ns ) ' 

(28) 



Y/i = S"(e ) denotes the second derivative of the self- 
energy operator over the energy argument, E"(e a ) = 
d 2 /{d 2 e)T 1 (e)\ £=e , V = T'(e a ) denotes the derivative 
of the vertex operator over the energy argument, and a' 
denotes the intermediate electron states with the energy 
£V = £a- The last term in Eq. (|2~7T) is added artifically, in 
order to make the whole expression for Ai?d.vr infrared 
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(IR) finite. The same term will be subtracted from the 
derivative contribution defined below, see Eq. (|29p . 

We note that all terms in Eq. (|27j) are UV finite, so that 
there is no need for any UV regularization. There are, 
however, IR divergences, which appear when the energy 
of the intermediate electron states in the electron prop- 
agators coincides with the energy of the reference state. 
The divergences cancel out in the sum of individual terms 
inEq. ([271). 



4- Derivative contribution 

Finally, the remaining contribution will be termed as 
the derivative term. It is given by 



AShht — 2 



/oo 



(aa'\I(u)\a'a) 
{-u + iO) 2 



x(V Zi 



(Sa ~ H)' 



Vhfs 



(29) 



The second term in the brackets is added artifically, in 
order to compensate the IR reference-state divergence 
present in the first term, making the total expression for 
AE^der IR finite. Note that this term is exactly the same 
as the one added to Eq. (|2"T1) but has the opposite sign. 

Finally, the total self-energy correction is given by the 
sum of the contributions discussed above, 



A£?se — A_Ep + AE VI hf s + AE V 



dcr ; 

(30) 



which are given by Eqs. ([11]), (JUJ), dHJ) , ([27]), and 
respectively. 



B. Self-energy: Calculation 

The general formulas reported so far represent contri- 
butions to the energy shift. We now have to separate 
out the nuclear degrees of freedom and convert the cor- 
rections to the energy into corrections to the shielding 
constant. In most cases, this is achieved simply by using 
the substitutions (fT3]) . The double- vertex contribution, 
however, requires an explicit angular-momentum algebra 
calculation for the separation of the nuclear variables. 

We will see that most of the corrections to the shield- 
ing constant can be regarded as generalizations of the 
corrections already discussed in the literature. So, our 
present calculation will be largely based on the previous 
investigations of the self-energy correction to the Lamb 
shift [H, Hq| ; to the hyperfine structure (3lT-l33l| , and to 
the g factor [32T4341 ] . The double- vertex correction of 
the kind similar to that in the present work appeared in 
the evaluation of the self-energy correction to the parity- 
nonconserving transitions in Refs. (35L l36j. However, in 
this work we develop a different scheme for the evaluation 
of the double- vertex contribution based on the analytical 
representation of the Dirac-Coulomb Green function. 



1. Perturbed- orbital contributions 

The matrix elements of the self-energy operator are 
diagonal in the relativistic angular momentum quantum 
number k and the momentum projection /i. Because of 
this, the angular reduction of the pcrturbed-orbital con- 
tribution is achieved by the same set of substitutions f| 13[) 
as for the leading-order magnetic shielding. The result- 
ing contribution to the shielding constant is conveniently 
represented as 

Aa-p = a{ai/ 2 \Y> R (s a )\5 {2) ai/ 2 ) 

+a(S$ s a. /2 \X R (e a )\Sgla 1/2 ) , (31) 

where the first-order perturbations of the reference-state 
wave function are given by 



1 



n i 



-■a °n 



(ni/ 2 |F h f s |ai/ 2 ) , (32) 



and 



(m/ 2 |Kee|ai/ 2 ) , (33) 



and IS^a) is the standard second-order perturbation [37j 
of the reference-state wave function induced by both in- 
teractions, Vhfs and Vzee- Note that only the diagonal 
in k part of the perturbed wave functions contributes 
to AtJp . The calculation of the non-diagonal matrix 
elements of the self-energy operator is performed by a 
straightforward generalization of the method developed 
in Ref. [30| for the first-order self-energy correction to the 
Lamb shift. 



2. Hfs-vertex contributions 

The hfs-vertex correction to the energy shift ([2"3"]l for 
the reference state with j a = 1/2 can be converted to the 
correction to the magnetic shielding by the substitution 
(fIBl. The result is 



Act vr , hfs = 7T~ 

Z7T 



dco 



ai/2n 2 \I(uj)\n 1 <$iccai/ 2 )("i|Vhf s |"2) 
[e a —ui- ue ni ){e a - u -ite„ 2 ) 



(Vhfs) 



lot 
2^ 



v-^ (ai/ 2 n\I(u))\n5illai/ 2 ) 

(e a -.-us n r ' (34) 



where the covariant regularization of the ultraviolet (UV) 
divergences is implicitly assumed. The right-hand-side of 
Eq. ((M]) differs from the vertex and reducible parts of the 
self-energy correction to the hyperfine structure only by 
the perturbed wave function |<5zoea) in place of one of the 
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reference-state wave functions \a). The main complica- 
tion brought by this difference is that the perturbed wave 
function contains components with different values of the 
relativistic angular quantum number k. So, for the ref- 
erence state with n a = — 1, the perturbed wave function 
has components with n = — 1 and k = 2, both of which 
contribute to the first term in Eq. denoted in the 

following as Aa VCI .Ms- The second (reducible) term con- 
tains only the k = K a component of the perturbed wave 
function and its calculation is done exactly as described 
in Ref. [3l[. Below, we present the generalization of for- 
mulas derived in Ref. [3l|, HH needed for the evaluation 

Of Acr V or,hfs- 

As explained in Ref. [31] , the covariant separation of 
the UV divergences is conveniently performed by dividing 
the vertex contribution into the zero- and many-potential 
parts, according to the number of interactions with the 
binding Coulomb field in the electron propagator, 



A(7 V 



= Act 



(0) 

ver.hfs 



T d+) 

ver.hfs ' 



(35) 



The zero-potential part is calculated in momentum space 
with the dimensional regularization of the UV diver- 
gences. The Fourier transform of the Vhts is done by 



(36) 



where q = p\ — p 2 is the transferred momentum. The 
contribution of the zero-potential hfs vertex part to the 



shielding constant is 



Act. 



(o) 



.hfs 



dpi dp 2 



X ^a 1/2 (Pl) 



(2tt) 3 (2tt) 3 
q x T R (px,p 2 ) 



^a 1/2 (p 2 ), (37) 



where p\ and p 2 are 4-vectors with the fixed time com- 
ponent pi = (e a ,Pi), P2 = (s a ,p 2 ), ip a and ip Sa are the 
reference-state and the perturbed wave functions, respec- 
tively, ip — ^^7° is the Dirac conjugation, and Tr is the 
renormalized one- loop vertex operator [28| . For evalu- 
ating the integrals over the angular variables, it is con- 
venient to use the following representation of the vertex 
operator sandwiched between the Dirac wave functions 



*Pa(Pi) r R (pi,p 2 )MP2) = 4- Kixl a „ a (Pi)vX-K b » b (jfo) 

+ ^2X- KaMa (Pl)0 ; X KbMtl (p2) 

+ (^3 Pi + K4p 2 )xi aUa {Pi)XK b v b (P2) 



KafJ>a 

(K5Pl+Tl6P2)X- Kafla (j>l)X-K b » t (j>2) , (38) 



where p = p/\p\, Xk^{p) are the spin-angular Dirac 
spinors [Hj], and the scalar functions TZi are given by 
Eqs. (A7)-(A12) of Ref. [3{|. Integration over the angu- 
lar variables yields (cf. Eq. (30) of Ref. (H) 



Act. 



(o) 

vcr.hfs 



48^ ^ ' 



Ja-h 



dpi r dp 2r 



Plr+P2r 



dq, 



P\rPlr 



|pir— P2r 



X I Tlx \pirKi(K a ,-Kg a ) -p 2r K[(K a , -KSa)} + U 2 \pi r Ki(-K a , K Sa ) - p 2r K[ (~K a , K Sa )} 
+ PlrP2r(7l 3 +K4) K 2 (n a ,K Sa ) + £>lrP2r C&5 + 7^6 ) K 2 (-K a , -«5o)} , 



(39) 



where pi r — \pi\, q r — \q\, nsa is the relativistic angu- 
lar quantum number of the perturbed wave function, 



+ 1/2 1 - 1/2, 



-2/3, x K=2 = -V2/3, 



and the basic angular integrals Ki(n, k') are defined and 
evaluated in Appendix fAl 



variables is 



Act. 



(1+) 

jr.bfs 



E 



ia 



dio 



K Sa = -l,2 

x X/ Y,j x j{ k i> k 2) Rj{u,an 2 ni8a) 



n 1 n 2 



(e - lj - ue ni ){e - ui - ue n2 ) 



subtraction , 



(40) 



The many-potential vertex contribution is free from 
UV divergences and thus can be calculated in coordinate 



(— 2 ) 

space. The result after the integration over the angular where R nin2 is the radial integral of the hfs type given 
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by Eq. (|T6|) . Xj is the angular coefficient, 



Xj(ki, k 2 ) 



]Sa ja J 



V2 

Kl — K 2 

~7T~ 



Ci(-K 2 , Hi) , 



(41) 



Ci(K a ,Kf,) is the reduced matrix element of the normal- 
ized spherical harmonics given by Eq. (CIO) of Ref. [28| . 
Rj is the relativistic generalization of the Slater radial 
integral given by Eqs. (C1)-(C9) of Ref. [2j|, and the 
subtraction in the last line of Eq. (PUfl) means that the 
contribution of the free propagators (already accounted 
for by the zero-potential term) needs to be subtracted. 



3. Zeeman-vertex contribution 

The Zeeman-vertex correction to the energy shift (|26|) 
for the reference state with j a — 1/2 can be converted to 
the correction to the magnetic shielding by the substitu- 
tion (IT3"1) . The result has the form analogous to that for 
the hfs-vertex contribution, 



Act v 



1 a 
2^ 



E 



{ai/ 2 n 2 \I{uj)\ni S^ayJ) (ni\V zest \n2 
(e a -oj- ue ni )(e a -ui- ue n2 ) 



IV ) — 



(ai /2 n\I(uj)\nd ( h 2ai/: 
(e a - ui -ue n ) 2 



duj 



(42) 



where the covariant regularization of the UV divergences 
is implicitly assumed. The above expression looks very 
similar to Eq. (f34|) and can be evaluated almost in the 
same way, except for the zero-potential vertex contribu- 
tion. The expression for the zero-potential Zeeman ver- 
tex contribution is different from the hfs case because 
the momentum representation of the interaction with a 
constant magnetic field involves a (5-function. The diag- 
onal matrix element of the Zeeman vertex operator was 
evaluated previously in Ref. [34| ; here we present the gen- 
eralization of the formulas required for the non-diagonal 
case. 

The Fourier transform of the interaction with the ex- 
ternal magnetic field is given by 



B x r 



-i(2ir) 3 B x V P 8 3 (p-f). 



(43) 



The contribution of the zero-potential vertex to the 
shielding constant is 



vcr.zec 



-7 CI 



dp dp' — 



(2tt) 3 



V p ~,6 3 (p~ 0) x f R (p,p') VwCpO , 
Jo' 2 

(44) 



where the gradient acts on the J-function only. This 
expression is transformed by integrating by parts and 
carrying out the integration with the J-function analyti- 
cally. The result after the angular integration consists of 
two parts (cf. Eqs. (27) and (36) of Ref. (U), 



vcr.zec 



where 

A (o) 

vcr.zec. 1 — _ /_j Jjk 

7T ' * 



A (0) A (0) 

vcr.zec, 1 vcr,zee,2 5 



(45) 



oo 2 



p dp 



A(p) 



o 

"•oa 

9a9Sa Al(K a , KSa) + fafSa Al(—K a , —Kg a ) 
P9afSa A 2 (K a , -KSa) ~ P faQSa A 2 (-K a , K Sa ) 

(46)" 



and 



Act. 



(o) 

ver,zee,2 



Of \ - 

\ ; 

Att £-^t ' 



Ala-h 



4TT ' ' — In »7T U 

KS a 



oo 2 



p dp 



X y>l{p) [g a fs a A2(K a ,-K Sa ) + fagSa A 2(-Ka,K8a)] 

+ ~ b l(p) [g a fSaA 3 (K a , -KSa) + fagSaA 3 (-K a ,K S a)] 
p 

+ h(p) [g a gSaAi{K ai K Sa ) + IafsaAi{-K a , ~K Sa )] 

+ h{p) [gagSaAi{K a , K Sa ) ~ fafSaA 4 (-K a , ~K Sa )] j , 

(47) 

where 5a„ = e a gsa+pfs a , fsa = £ a fsa+pgs a , g a = £ a g a + 
Pfa, fa = £afa+pg a , g' = dg(p) / 'dp and /' = df(p)/dp, 
the scalar functions A(p), bi(p) are given by Eqs. (24), 
(30)-(32) of Ref. [34{ and Ai(ni, k 2 ) are the basic angular 
integrals defined and evaluated in Appendix [X] 



Note that the expression for Act 



(o) 

vcr.zec, 2 



is non- 



symmetric with respect to a and 5a, but the result 
does not change when the two wave functions are in- 
terchanged. 



4- Double-vertex contribution 

The double- vertex correction is defined by Eq. (j27|) . 
All parts of it are UV finite and thus can be evaluated in 
coordinate space. We denote the individual terms in the 
right-hand-side of Eq. (|27f by SiE and consider each of 
them separately, 



AE, 



d.vr 



SiE. 



(48) 



The first term is 



5iE = 2 (A) = 2 



2tt 



duj 



{anz\I{u})\n\a){nx\ V zee \n 2 ){n 2 \ V his \ n 3 ) 



n^n 3 ~ u - ue ni ){e - u - ue n2 )(e - u - ue n3 ) ' 

(49) 



After separating the nuclear variables and integrating 
over the angles, the contribution to the magnetic shield- 
ing is (for j a = 1/2) 



• 2 poo 

J-oo ni n 2 n 3 J 

(1) d("2) 



Rj(uj, an 3 n\a) 



(e-w- ue ni )(e - uj - ue n2 )(e - uj - ue n3 ) ' 

(50) 



where 



X d j VCI ( Kl ,K 2 , K3 ) 



f—iy+3a-ja 






jn 


3 a 1 




3a 1 




3i 


32 / 




J2 J 



x £ (2j„ + 1) 

j„=l/2,3/2 
«1 + K 2 K 2 + K 3 

x - — Oi(-«2,*siJ j= — Gi(-k 3 ,k 2 ). (51) 



V3 



V3 



Note that the expression for 5\a is IR divergent when 
any two or all three of the intermediate states have the 
same energy as the reference state: s ni = e„ 2 = e a , 

£n-2 £-a-> ^-ri3 — £-a-> — ^ri2 — ^ri3 — 

The IR divergence cancels when all parts of Eq. (|27[) are 
added together. 

The contribution of the second term to the shielding is 

S 2 a = f (Kee) (V his ) <E") = ~ Rg> RL 2) 

ia 2 f°° , 



x du y 

2t ./_«, ^ ^ 2j a + 1 



,K2 ^Kl 



Rj(ui, an^nxa) N nin2 N n2n3 



(e — u — ue ni )(e — U3 — ue ri2 )(e - uj - ue„ 3 ) ' 

(52) 



where N a b is the normalization integral, 



N ab = / dxx {g a gb + fafb) ■ 
Jo 



(53) 



The expression for 5 2 a is IR divergent when s ni = e n2 

£-ri3 — 



The third term is given by 



J2j Xj(ki, k 2 ) 6 K2 ^ 3 Rj(uj, an 3 nia) R 

111 '12 



(e - uj - ue ni )(e - uj - ue 1l2 ){e - uj - us na ) 
J2j Xj(k 2 , Kg) 6 Kl , K2 Rj(uj, an 3 nia) N niTl2 R 



(-2) 
ri2 "3 



(e - ui - ue ni )(e - uj - ue n2 )(e - uj - ue n3 ) 



where 



Xj{k\, k 2 ) 



1 

71 



h h 1 

ja ja 'J 



-Kl 



«2 



V3 



(54) 

Cl(-K 2 ,Kl) . 

(55) 



The fourth term 8^a is obtained from 5 3 a by an obvious 
substitution O i?<~ 2 ). 
The fifth term is given by 



2^ 



aw -, —tt. . (56) 

{-uj + tO) 2 v ; 



Finally, the total double-vertex contribution is given 
by the sum of the five terms discussed above, 



AfJd.vr = ^ 6j<J • 



(57) 



i=l 



Despite the fact that the individual terms 8i<7 are IR di- 
vergent, the sum of them is finite and can be evaluated 
without any explicit regularization, provided that (i) the 
integration over the frequency of the virtual photon in 
the self-energy loop ui is performed after all five terms are 
added together and (ii) the contour of the uj integration is 
suitably chosen. One can show that if the uj integration is 
performed along the contour consisting of the low-energy 
and high-energy parts, as e.g., in Refs. [13,|3l|, the inte- 
grand becomes regular at ui — > and, therefore, can be 
directly evaluated numerically. 

The numerical evaluation of the double- vertex correc- 
tion Acd.vr was the most time consuming part of our 
calculation, due to a large number of the partial waves 
involved and a four-dimensional radial integration in 
Eq. (l50l) . The radial integration was carried out with 
help of the numerical approach developed in our calcula- 
tion of the two-loop self-energy and described in detail in 
Ref. [2!| . Because of numerical cancellations between the 
five terms in Eq. (|57|) . especially in the region of small 
values of w, we took care to treat all five terms exactly 
in the same way. In particular, the normalization inte- 
grals N a b in Eqs. (j52l and (f5~4"]l were evaluated numer- 
ically, in order to be consistent with the evaluation of 
Eq. (|50l) . whereas the corresponding contributions could 
have been evaluated more easily by taking the derivative 
of the Dirac-Coulomb Green function. 
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C. Vacuum polarization correction 

Vacuum polarization corrections to the magnetic 
shielding calculated in the present work are shown in 
Fig.[U The diagrams in Fig.[3Ja)-(c) come from the inser- 
tion of the Uehling potential into the electron line of the 
leading-order magnetic shielding, whereas the diagram 
in Fig. EJd) respresents the Uehling potential insertion 
into the hfs interaction. We note that the diagram with 
the vacuum polarization insertion into the Zeeman inter- 
action vanishes in the Uehling-potential approximation. 
In our treatment, we neglect contributions with addi- 
tional Coulomb interactions in the vacuum-polarization 
loop (which correspond to the Wichmann-Kroll part of 
the one-loop vacuum polarization) and the additional di- 
agram of the Wichmann-Kroll type with both hfs and 
Zeeman interactions attached to the vacuum-polarization 
loop. We expect that the part accounted for in the 
present work yields the dominant contribution to the vac- 
uum polarization correction. 

The contribution of the diagrams in Fig. [5] (a)-(c) is 
an analogue of the perturbed-orbital self-energy contri- 
bution and is given by a similar expression, 

Acrvp,po = a (ai /2 |Vueh|<5 (2) ai/ 2 ) 

+a (4feOy 3 | Vbehl^ai/2) 
where Vijoh is the Uehling potential, 
2a 2 Z 



(58) 



U Ueh (r) 



dr r p(r ) 



x [K (2m\r - r'\) - K (2m\r + r'\)} , (59) 



and 



K (x) = 



dtt 



V + w)^' (60) 



and the nuclear-charge density p is normalized by the 
condition J dr p(r) = 1. We note that this contribu- 
tion can also be calculated by incorporating the Uehling 
potential into the Dirac equation and re-evaluating the 
leading-order magnetic shielding, in this way account- 
ing for the Uehling potential to all orders. We performed 
calculations in both ways, which ensured that the pertur- 
bations of the reference-state wave function \S^a) and 
|j' 2 'a) are computed correctly. 

The contribution of the diagram in Fig. [U(d) to the 
magnetic shielding can be expressed as 



Acryp,: 



where 



x (ai/2|V r Z oo|ni/ 2 )(ni /2 |U V p, mag |ai/ 2 ) , (61) 



Uvp imag (r) 



Vhft(r) 
1 + 



2a 

3tt ji 
1 

2*2 



dt 



Vt 2 



t 2 



(1 + 2mrt) e" 



-2rnrt 



(62) 



is the hfs interaction modified by the vacuum- 
polarization insertion. 



IV. QED CORRECTION FOR SMALL 
NUCLEAR CHARGES 

In the previous section, we calculated the QED correc- 
tions to the magnetic shielding without any expansion in 
the nuclear binding strength parameter Za. Now we turn 
to the evaluation of this corrections within the expansion 
in this parameter. We will derive the complete expres- 
sion for the leading term of the Za expansion, which 
enters in the relative order 5a /a ~ a{Za) 2 . The re- 
sults obtained in this section will be applicable for light 
hydrogen-like ions, where no all-order calculations were 
possible because of large numerical cancellation. They 
will also provide an important cross-check with the all- 
order calculations described in the previous section. 

For the derivation, it is convenient to use the for- 
malism of the non-relativistic quantum electrodynamics 
(NRQED). In this approach, all QED effects are calcu- 
lated by an expansion in powers of a and Za and repre- 
sented as expectation values of various effective operators 
on the nonrelativistic reference-state wave function. Let 
us start with the effective Hamiltonian -//nrqed, which 
includes leading one-loop radiative corrections: 



-Hnrqed = he A -- . , 

2 m V 4 m z 



3 K 
4 m 2 r 

Q-M ) B 2 



12 m 

„2 



{7?, VxJj} 



4 TO 3 



(63) 



where { . , .} denotes anticommutator, 7? = p — e A. The 
QED effects in the above Hamiltonian are parameter- 
ized in terms of the constants k, Te, and a m, which are 
interpreted as the electron magnetic moment anomaly, 
the charge radius and the magnetic polarizability, respec- 
tively, 



OtM = 



a 

2^' 
3 K 
z m z 
4a 
3 7T m 3 





' m 


2 ° ( 


In — - 


7T m z \ 


v 2e 




13\ 




+ 2l) 


V 2e 





-) 

24/ 



(64) 
(65) 

(66) 



where e is the photon momentum cutoff. In the non-QED 
limit, all QED constants vanish, r 2 E — olm = k = 0, 
and the effective Hamiltonian -//nrqed turns into the 
Schrodinger-Pauli Hamiltonian. 

The QED constants r 2 E and om depend explicitly on 
the the photon momentum cutoff parameter e. This de- 
pendence cancels out with contributions coming from the 
emission and absorption of virtual photons with the en- 
ergy smaller than e. The complete expression for any 
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(a) ^ 



(b) | 



(c) ^ 




FIG. 2: Vacuum-polarization corrections to the nuclear magnetic shielding calculated in the present work. The double line 
represents electron propagating in the binding nuclear field, the single line represents the free-electron propagator, the dashed 
line terminated by a cross represents interaction with the binding Coulomb field, the wave line terminated by a triangle 
represents the dipole hyperfine interaction with the nucleus, and the wave line terminated by a cross represents the interaction 
with the external magnetic field. 



physical quantity, e.g., the Lamb shift and the shielding 
constant, does not depend on the artificial photon cutoff 
parameter. 

The above formula for r 2 E is obtained from the well 
known expressions for the one-loop radiative corrections 
to electromagnetic formfactors F\ and F2 [1(3]. The for- 
mula for olm has not appeared in the literature previ- 
ously. We have derived it by a method similar to that 
used in Ref. [4l[ for the electric polarizability, denoted 
by x m that work. 

We will start with rederiving the known nonrelativis- 
tic expression for the shielding constant within our ap- 
proach. It comes from the A 2 term in the electron kinetic 
energy in Eq. (|63|) . The electromagnetic potential A is 
the sum of the external magnetic potential Ae, 



A E = - B x f 
2 



(67) 



and the potential induced by the nuclear magnetic mo- 
ment, 



7 1 - ** 

M = -;— i* x -3 • 

4 7r r° 



(68) 



The corresponding energy shift is 
SE 

2m 1 

' B x r 



P 2 -*2 P 2 -* -* 

— (A ) = -(A E -Aj) 
2m m 



a 
2 m 



r 

V x ^3 



(69) 



where matrix elements are calculated with the nonrela- 
tivistic wave function. The shielding constant a is ob- 
tained from the energy shift 5E by 5E = a fx ■ B. For 
the ground (L = 0) hydrogenic state, the nonrelativistic 
results is 



a I 
3 m \ r 



Before considering the QED effects to the shielding 
constant, it is convenient to first recalculate the leading 



QED correction to the energy levels. The total contribu- 
tion is split into two parts induced by the virtual photons 
of low (L) and high (H) energy, 



SE Lamh = SE L + 5 Eh 



(71) 



The high-energy part is the expectation value of the V • E 
term in the effective Hamiltonian in Eq. (|63l) 



a (Za) 4 
it n 3 



4 , m 10\ r 

— In 1 O/o 

3 2e 9 ' 



(72) 



where E = -VA° and A = -Z e/(4irr). The vacuum 
polarization can be incorporated in the above expression 
by adding — 2a/(5 7rm 2 ) to r% in Eq. (IB31) . The low- 
energy part SEl is induced by the emission and the ab- 
sorption of the virtual photons of low (fc < e) energy, 



SEr 



d 3 k 



2a 



3 7r \ m 



P 



o (2n) 3 2k 
(H-E) jln 



g*3 



P' 



1 



V 



m E-H-k 



2 c 



m (Z a) 



In 



2 (H-E) 
m (Z a) 2 



(73) 



The terms containing lne cancel out in the sum of the 
low- and high-energy parts, as expected. The total 
leading-order Lamb shift contribution is 



a (Zaf [4, 

O-C/Lamb = 3 — \ 77 in (Z a) 

7T il 3 



10 4 ^ A 4 1 fcC n 
o/n In kn (n, t 

9 15/ 3 1 ' ' 



where lnfco(n, I) is the Bethe logarithm, 

3 *2(H-E) 



(74) 



(70) hifco(M) 



2m 3 (Za) 4 



p(H - E) In 



m (Z a) 2 



(75) 

We now turn to the calculation of the QED correction 
to the magnetic shielding, which is performed similarly 
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to that for the Lamb shift. The total contribution is split 
into the low and the high-energy parts: 



where 
SE L = 

and 
SE H -- 



8E = SE L + SE H 
d 3 k 



(76) 



o (2 7r) 3 2fc 



S IJ 



r E V-E 



1 



k l k 3 



.4' 



7T* 1 TP? 

to E — H — k to 



(77) 



12 TO 



(E — H)' 2 m 

Pi) {-,VxB} 
4 m l I k ' 



e -2 
y a« B 



The high-energy part can be conveniently rewritten in 
the form 



1 



1 



7 4?r5 d (r) )fi-B 



2 7r a 



3 K 



(I £ -Vx Bi) 



3 m \ ^ 4 to 2 

-4:Traa M (B ■ Bj) . (79) 

In the low-energy part, we separate out lne and then 
perform an expansion in the magnetic fields, 

5E L = SE LA + SE LB , 



3 7T \ TO TO(Ztt) z TO 

^ = |^(tf-i^)ln-^. 

3 7T \ TO TO/ TO (Z a) 

Using the identity 



(80) 
(81) 

(82) 



2 - (B - B) - 

. TO TO 



7T 

m 



H - E, — 

TO 



4 7rZa(5 3 (r) + -t^Itt. V x B j + 

2 m J I J to 



(83) 



we transform SElb to the form 



SE LR = — In 



3vr 



2r 



to (Z a) 2 
1 



- At: Z aS (r) j fl ■ B (84) 



3r {E-H] 
Air a (A E ■ V x B/) + 167ra (B • Bj) 



The artificial parameter In e cancels out in the sum 
<5i?LB + SEh separately for each type of the matrix ele- 
ments, 

2 a 2 

6E LB + 5E H = - — Za 

9 7T 

5 1 



6 



fi- B (Air8 6 {r) 



3° 
1G 



31 1 



"ln(^«)- 2 -|" 



E-H)' r 
(Ae • V x Bj) 



(B • Bj). 



(85) 



Using the following results for the matrix elements with 
nS states 



47r<5 3 (r) 



E-H) 1 r 
(Ie'Vx Bj) 

(B • Bj) 



6(Za) 2 

3- M ' B, 

7r 

o J— M'B, 

3 7r n J 



(86) 
(87) 



we obtain 



(78) 6E LB + SE H = — 



a 2 (Za) 



9 7r n 3 



, , x 2 421 3 

n(ZaU 2 + - 

V ' 96 5 



M-B. 
(89) 

The calculation of the remaining low-energy contribu- 
tion 8Ela is slightly more complicated. We first return to 
the integral form of 8 El, derive an expression for 8 El a, 
and then drop all terms with lne, as they are already 
accounted for by 8Elb- The integral form of 8El is 



SE L = ?H / /,■ ,//,- / _ 



3tt 7o 



to E — H — k to 



(90) 



It can be rewritten, using the identity 7? = — i to [r , H — 
E] , in the form 



8Er 



2a 
3^ 



k 3 dk(r 



E-H -k 



(91) 



All terms with positive powers of e are discarded, since 
one assumes that the limit e — > is taken after the ex- 
pansion in a is done. We shall now expand the integrand 
in Eq. (|9"Tj) in the magnetic fields. The Hamiltonian H 
is 



ry -* f> -> -t €- -> 

H = Ho + — fl ■ B — —— L ■ B — ji'Mi (92) 

3 r 2 TO 4 7T 77Z T 



where Bo = p 2 /(2 m) — Za/r. The first term with B in 
Eq. jS2} can be absorbed into Z' = Z - fl ■ B/3. So, the 
correction due to this term is 



SE L i 



{Z'af \ In 



4a 
3?r 

16 a 2 (^a) 3 
9 7r n 3 



2c 



to (Z' a) 2 



In fcg 



(93) 



1 



fl ■ B < In fco + — — In 



2r 



to (Z a) 2 



The correction due to the two last terms in Eq. (|92j) is 
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SE L2 = 



2 a 2 

37TTO 2 



dkk S 



1 



2 a 2 
9tt 



p- B 



dk k 



E-H-k 



LB 



1 



dk\ E-H 



E-H-k r 3 

1 1 
k r 3 E-H-k 



E-H-k 



(94) 



Integrating by parts and using the results and the nota- 
tion of In &3 from Ref. [l2| 



dkk 2 { r 
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E-H- kr 3 E-H-k 
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2e 
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m [Z a) 2 


-hifcsj 







(95) 



we obtain 



6E L o = 



8 a 2 {Zaf 



3 7r m 2 n 3 





2e 




{'•• 


m (Z a) 2 


-m*s-|} 







(96) 

The sum of these two contributions SEla = SEh+SEl2, 
after dropping the lne terms, is 

5£LA = ^^^M-5(2 1nA;o-3 1n/c 3 ). (97) 

Finally, the total correction to energy is 

SE = 5E LA + 5E LB + 5E H . (98) 

Expressing the energy shift in terms of the shielding con- 
stant (SE — 8a fl ■ B), we obtain the complete result 
for the leading QED correction to the nuclear magnetic 
shielding in hydrogen-like ions, valid for the nS states, 



8 a 2 {Za) 3 

6(7 = 3~~ 

9 7r n J 



(99) 



f 421 3 

: lln[(Za)- 2 ] + 2 In k {n) - 3 In k 3 {n) - — + - 

(100) 



The term of 3/5 in the brackets is the contribution of 
the vacuum polarization. The numerical results for the 
Bethe logarithm In ko and the 1 /r 3 Bethe- logarithm-type 
correction lnfc3 [l2[ for the Is state are 



lnfco(ls) = 2.984128 556, 
lnfc 3 (ls) = 3.272 806 545. 



(101) 
(102) 



We note that the numerical value of the constant term 
in Eq. (|100p . —7.635 58, is comparable in magnitude but 
of the opposite sign to the logarithmic term at Z = 1, 
lna~ 2 = 9.840 49. This entails a significant numerical 
cancellation between these two terms for hydrogen and 
light hydrogen-like ions. As a result, the total QED cor- 
rection turns out to be much smaller in magnitude than 
could be anticipated from the leading logarithm alone. 



V. OTHER CORRECTIONS 
A. Bohr-Weisskopf correction 

We now turn to the effect induced by the spatial dis- 
tribution of the nuclear magnetic moment, also known 
as the Bohr-Weisskopf (BW) correction. Our treatment 
of the BW effect is based on the effective single-particle 
(SP) model of the nuclear magnetic moment. Within this 
model, the magnetic moment is assumed to be induced 
by the odd nucleon (proton, when Z and A are odd and 
neutron, when Z is even and A is odd) with an effec- 
tive g-factor, which is adjusted to yield the experimental 
value of the nuclear magnetic moment. The treatment 
of the magnetization distribution effect on hfs within the 
SP model was originally developed in Refs. [42|, HH and 
later in Ref. (44|. Our present treatment closely follows 
the procedure described in Refs. (4#t46l|. 

The wave function of the odd nucleon is assumed to 
satisfy the Schrodinger equation with the central poten- 
tial of Woods-Saxon form and the spin-orbital term in- 
cluded (see, e.g., Ref. [47j]) 

V(r) = -V Q F(r) + — so (r)T- a + V c {r) , (103) 



where 



4>so(r) 



V so dF{r) 



4m p r dr 



i 



exp 



r-R 



(104) 



(105) 



and Vc is the Coulomb part of the interaction (absent 
for neutron) , with the uniform distribution of the charge 
(Z — 1) over the nuclear sphere. The parameters Vq, V so , 
R, and a were taken from Refs. (47ll48|. 

The nuclear magnetic moment can be evaluated within 
the SP model to yield ^ 



JL 

fJ-N 



1 

2 93 



I- 1 - 
2 



= < 



2(7+1) 



9s 



21+1 
4(/ + l) 



7(27 + 3) 



2(7+1) 



9l 

for I - 

27+1 
4(7 + 1) 

for I - 



9l ■ 



1 

~ 2 ' 
(106) 
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where 7 and L are the total and the orbital angular mo- 
mentum of the nucleus, respectively, g^ is the g factor 
associated with the orbital motion of the nucleon (<?l = 
1 for proton and g^ = for neutron) and gs is the effec- 
tive nucleon g factor, determined by the condition that 
Eq. (|106|) yields the experimental value of the magnetic 



moment. 

It was demonstrated in Ref. [45[ that, within the SP 
model, the BW effect can be accounted for by adding a 
multiplicative magnetization-distribution function F(r) 
to the standard point-dipole hfs interaction. The distri- 
bution function is given by [4(3] 



F(r) 



UN 

H Jo 

| UN 

n 



dr' r' 2 \u(r')\ 2 



9S 



I 



1 27 + 1 



r 2 <j> so (r) gh 



dr' r' 2 \u{r')\ 2 



r.3 r 



2 4(7+1) 
27-1 / 1 27 + 1 



8(7 + 1) 



gs 



i - 



2 4(7+1) 



r (/>so(r) gh 



for 7 = L + 1/2 and 

F(r) = ^ r dr'r' 2 \u{r 
M Jo 



'M2 



2(7+1) 



gs 



1(21 + 3) 27+1 2 

^(7TiT-4(7TT) r ^ (r) l 5L 



Mat / , / ,2, , ,s |2 r 



27 + 3 / 7(27 + 3) 27+1 



8(7+1) 55 ' V 2(7+1) 4(7+1) 



r <l>so{r) gh 



(107) 



(108) 



for 7 = L — 1/2. In the above formulas, u(r) is the wave 
function of the odd nucleon. It can easily be seen that 
F(r) = 1 outside the nucleus. 



B. Recoil and quadrupole corrections 

The recoil correction to the magnetic shielding was ob- 
tained in Ref. [18| in the nonrelativistic approximation, 



6 (Jr. 



aZa m 
3 M 



9N-1 
9n 



where M is the nuclear mass and 

M (i 



gN 



Zrrip hn 7 



(109) 



(110) 



The electric-quadrupole correction to the magnetic 
shielding is 



3Q 



27(27- 1) (m/m p )g I 



(111) 



where Q is the quadrupole moment of the nucleus and 
SgQ is the quadrupole correction to the g factor calcu- 
lated in Ref. fijl, 



Sg Q = a (Zaf 



= a(Zaf 



12[35 + 20 7 - 32(Za) 2 ] 

135 7(l + 7) 2 [15- 16(Za) 2 

11 43 . , 2 

1 (Za 2 + . 

135 405 v 1 



(112) 



VI. RESULTS AND DISCUSSION 

Numerical results for the self-energy correction to the 
nuclear magnetic shielding can be conveniently parame- 
terized in terms of the dimensionless function Dg^(Za) 
defined as 



AdsE = a 2 (Za) 3 D SE (Za) 



(113) 



Our numerical all-order (in Za) results for the self-energy 
correction to the magnetic shielding are summarized in 
Table fl] We note significant numerical cancellation be- 
tween the individual contributions, which is particularly 
strong for low values of Z. The fact that the result- 
ing sum is a smooth function of Z and demonstrates the 
expected Z scaling serves as a consistency check of our 
calculations. Besides of numerical cancellation, the ad- 
ditional complication arising in the \ow-Z region is that 
the convergence of the partial-wave expansions becomes 
slower when Z decreases. Because of these two compli- 
cations, the accuracy of our results worsens for smaller 
values of Z and no all-order results were obtained for 
Z < 10. 

The all-order numerical results can be compared with 
the Za-expansion results obtained in Sec. IIVI As follows 
from Eq. (|100[) . the Za expansion of the function 7?se 
for the Is state reads 



D SE (Za) = 



9tt 



\n(Za)' 2 - 8.235 579 + 0(Za) 



(114) 



where O(Za) denotes the higher-order terms. In Fig. |31 
the numerical all-order results for the function 7?se are 
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FIG. 3: (Color online) Self-energy correction to the nuclear 
magnetic shielding. 




Z 

FIG. 4: (Color online) Vacuum-polarization correction to the 
nuclear magnetic shielding. 

plotted together with the contribution of the leading log- 
arithmic term in Eq. (II 14[) (dashed line, red) and the 
contribution of both terms in Eq. (|114[) (dashed-dot line, 
blue) . We observe that the leading logarithm alone gives 
a large contribution that disagrees strongly with the all- 
order results. However, when the constant term is added, 
the total Za-expansion contribution shrinks significantly 
and even changes its sign for Z > 3. Only after the con- 
stant term is accounted for, we observe reasonable agree- 
ment between the all-order and Za expansion results. 

The vacuum-polarization correction to the nuclear 
magnetic shielding is parameterized as 

(Za) 3 Dyp(Za) . (115) 

Our numerical results for the vacuum-polarization cor- 



FIG. 5: (Color online) Individual contributions to the nuclear 
shielding. "NR" is the nonrelativistic contribution, "REL" 
is the relativistic point-nucleus contribution, "FNS" is the 
finite nuclear size correction, "QED" is the QED correction, 
"BW" is the Bohr-Weisskopf correction, "REC" is the recoil 
correction, and "QUAD" is the electric quadrupole correction. 
Note that the QED correction changes its sign between Z = 4 
and 5. 

rection are presented in Table |TIJ The calculation was 
performed for the extended nucleus and the Uehling po- 
tential was included to all orders. Note that our present 
treatment is not complete, as the Wichmann-Kroll part 
of the correction is still missing. We estimate the un- 
certainty due to omitted terms to be within 30% of the 
calculated contribution. 

The Za expansion of the function -Dyp reads 

D VP (Za) = -^^5 l , + O{Za) ) (116) 

where the Kronecker symbol Sifi indicates that the cor- 
rection vanishes for the reference states with I > 0. Com- 
parison of the all-order numerical results with the lead- 
ing term of the Za expansion is given in Fig. 01 It is 
remarkable that the all-order results grow fast when Z 
is increased and eventually become more than an order 
of magnitude larger than the leading-order result. In the 
high-Z region, the self-energy and vacuum-polarization 
corrections are (as usual) of the opposite sign, the self 
energy being about twice larger than the vacuum polar- 
ization. In the \ow-Z region, however, the self-energy 
correction changes its sign and the total QED contribu- 
tion becomes positive for very light ions. 

The summary of our calculations of the nuclear mag- 
netic shielding constant a for several hydrogen-like ions 
is given in Table IIIII The first line of the table presents 
results for the leading-order nuclear shielding, including 
the finite nuclear size effect. The leading-order contri- 
bution was calculated using the Fermi model for the nu- 
clear charge distribution and the point-dipole approxi- 
mation for the interaction with the nuclear magnetic mo- 
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ment. The nuclear charge radii were taken from Ref. [49j. 
The results arc in good agreement with those reported in 
Ref. [17j. The QED correction presented in the second 
line of Table IIIII is the sum of the all-order results for the 
self energy and the vacuum polarization. Its error comes 
from the numerical uncertainty of the self-energy part 
and the estimate of uncalculated vacuum-polarization di- 
agrams. Since there were no all-order calculations per- 
formed for oxygen, we used an extrapolation of our all- 
order results (taking into account the derived values of 
the Za expansion coefficients). 

The Bohr-Weisskopf correction, presented in the third 
line of Table IIII| was calculated by reevaluating the 
leading-order contribution with the point-dipole hfs in- 
teraction modified by the extended-distribution function 
F(r) given by Eqs. (|107|) and (|108|) , Because the effective 
single-particle model of the nuclear magnetic moment is 
(of course) a rather crude approximation, we estimate the 
uncertainty of the Bohr-Weisskopf correction to be 30%, 
which is consistent with previous estimates of the uncer- 
tainty of this effect [45j . This uncertainty includes also 
the error due to the nuclear polarization effect, which is 
not considered in the present work. 

The quadrupole and the recoil corrections, given in 
the fourth and fifth lines of Table IIIII respectively, were 
evaluated according to Eqs. (|111[) and (|109|) . The er- 
ror of the quadrupole correction comes from the nu- 
clear quadrupole moments. The magnetic dipole and 
electr ic q uadrupole moments of nuclei were taken from 
Refs. [M[5l|. 

The Z dependence of individual contributions to the 
magnetic shielding constant a is shown in Fig. [5] The 
leading-order contribution is separated into three parts, 
the point- nucleus nonrelativistic part (solid line, black), 
the point-nucleus relativistic part (dashed line, green), 
and the finite nuclear size correction (dashed-dotted line, 
blue). We observe that the finite nuclear size correction, 
as well as the other effects calculated in this work, be- 
come increasingly important in the region of large nu- 
clear charges. The only exception is the nuclear recoil 
correction. It practically does not depend on the nuclear 
charge number Z, since the linear Z scaling in Eq. (| 109[) 
is compensated by the increase of the nuclear mass M 
with Z . As a consequence, the recoil effect is completely 
negligible for high- and medium- Z ions, but turns into 
one of the dominant corrections for Z < 10. 

For most ions, the uncertainty of the theoretical pre- 
diction is roughly 30% of the Bohr-Weisskopf effect; it 
can be immediately estimated from Fig. [SJ For very light 
ions, however, there is an additional uncertainty due to 
the unknown relativistic recoil effect. Its contribution 
can be estimated by multiplying the nonrelativistic re- 
coil correction plotted on Fig. [5] by the factor of (Za) 2 . 
We observe that for very light ions, the relativistic recoil 
is the dominant source of error in theoretical predictions. 



VII. CONCLUSION 

In this work we have performed an ab initio calcula- 
tion of the nuclear magnetic shielding in hydrogen-like 
ions with inclusion of relativistic, nuclear, and QED ef- 
fects. The uncertainty of our theoretical predictions for 
the nuclear magnetic shielding constant defines [accord- 
ing to Eq. ([5])] the precision to which the nuclear mag- 
netic dipole moments can be determined from experi- 
ments on the ^-factors of hydrogen-like ions. It can be 
concluded from Table ITO1 and Fig.[S]that the present the- 
ory permits determination of nuclear magnetic moments 
with fractional accuracy ranging from 10 -9 in the case 
of 17 O 7 +tol0- 5 for 209 Bi 82 +. 

For most hydrogen-like ions, the dominant source of 
error in the theoretical predictions is the Bohr-Weisskopf 
effect. Since this effect cannot be accurately calculated 
at present, we conclude that the theory of the nuclear 
magnetic shielding has reached the point where the un- 
certainty due to nuclear-structure effects impedes further 
progress. For very light ions, however, the dominant the- 
oretical error comes from the unknown relativistic recoil 
effect, whose calculation might be a subject of future 
work. 
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Appendix A: Angular integrals 

Throughout this work, we repeatedly used the follow- 
ing result for the matrix element of the hfs and Zeeman 
interaction 

M^^M =(-1)*+"" 

(Al) 

where r = f/\r\, „, is the Clebsch-Gordan coef- 

ficient, and 

P<*>(ni,n 2 ) = K2 dC-zca, *i) R^ n2 ■ (A2) 

The radial integral is defined by Eq. (fTB) and the 
angular coefficient Ci(«; a , «(,) is the reduced matrix el- 
ement of the normalized spherical harmonics, see, e.g., 
Eq. (CIO) of Ref. [2§j |. An important particular case is 
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TABLE I: Individual contributions to the self-energy correction to the nuclear magnetic shielding, for the point nucleus, in 
terms of the function D$e defined by Eq. (|113p . 



z 


po 


vr,hfs 


vr,zee 


d.ver 


der 


total 


10 


-9.584 


6.095 


9.180 


-49.722 


43.523 


-0.508(100) 


14 


-4.741 


1.711 


4.321 


-21.048 


19.047 


-0.710(15) 


16 


-3.571 


0.705 


3.140 


-14.718 


13.655 


-0.789 (9) 


20 


-2.217 


-0.405 


1.761 


-7.848 


7.782 


-0.927 (4) 


26 


-1.277 


-1.115 


0.772 


-3.473 


3.983 


-1.110 (2) 


32 


-0.858 


-1.407 


0.292 


-1.644 


2.333 


-1.283(1) 










U. JOu 


l.OU 


1 ^1 Q !~\ \ 


45 


-0.575 


-1.643 


-0.171 


-0.267 


0.975 


-1.681 


54 


-0.595 


-1.744 


-0.335 


0.023 


0.622 


-2.029 


60 


-0.672 


-1.825 


-0.423 


0.112 


0.487 


-2.321 


70 


-0.929 


-2.028 


-0.571 


0.175 


0.355 


-2.999 


82 


-1.612 


-2.508 


-0.812 


0.191 


0.284 


-4.457(2) 


83 


-1.701 


-2.569 


-0.839 


0.192 


0.281 


-4.636(1) 


92 


-3.011 


-3.400 


-1.174 


0.198 


0.280 


-7.107(2) 



TABLE II: Vacuum-polarization correction to the magnetic 
shielding, for the extended nucleus, in terms of the function 
Dy P defined by Eq. (|115|) . 



z 


po 


mag 


total 


10 


0.118 


0.110 


0.228 


14 


0.135 


0.121 


0.256 


16 


0.144 


0.126 


0.271 


20 


0.164 


0.137 


0.302 


26 


0.200 


0.155 


0.355 


32 


0.242 


0.175 


0.417 


40 


0.314 


0.206 


0.520 


45 


0.369 


0.228 


0.597 


54 


0.500 


0.275 


0.775 


60 


0.618 


0.315 


0.933 


70 


0.891 


0.398 


1.289 


82 


1.449 


0.546 


1.996 


83 


1.512 


0.562 


2.074 


92 


2.227 


0.727 


2.954 



K\ = k-2 = — 1 and fii = 
above formulas reduce to 



/i2 = V 2 ; m which case the 



-iy 2 | 



Ir x a) . 



li/ 2 ) 



-3^6 



(-a) 



(A3) 



The basic angular integrals Ki needed for the evalua- 
tion of the hfs vertex contribution are defined by 



3i ' 
4tt 



dpx dp 2 F(£) x Ka i /2 (Pi) [Pi x o]o X^iA) = / d€F(£) K x (K a , K b ) , (A4a) 



4tt 

47T 



-1 
1 



|^ f dfx d£ 2 F(0 xl aV M) $2 x a] X.M&) = / <% F (t) #i(«a, K b ) , (A4b) 



dpi dp 2 F(£)xi al/2 (pi) [pi x P2]qX^/2(P2) = / d^F(^)K 2 (Ka,K b ) , (A4c) 



where F(£) is an arbitrary function and ^ — Pi ■ P2- The integrals over all angles except for £ in the above equations 
are evaluated analytically, as described in Appendix of Ref. [33[ . The results relevant for this work are 

Ki(L2) = j- 2 , (A5) 



Ki(- 


-1,1) = 


-6 


^(1,-1) 


= 1, 


#i(- 


-1,-2) = 


K[(- 


-1,1) = 


-1, 


^(1,-1) 


= £, 


K[(- 


-1,-2) = 


K 2 {- 


-1,-1) 


= 0, 


K 2 (l, 1) = 


'2 ' 


K 2 (- 


-1,2) = 0, 
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TABLE III: Individual contributions to the shielding constant a x 10 for selected hydrogen-like ions. 



^O^ 43 Ca 19 + 73 Ge 31 + 131 Xc 53+ 209 Bi 82 + 

Leading 143.3127 375.960 657.93 1461.6 4112 

QED -0.0026(2) -0.103(15) -0.59(8) -4.1 (0.8) -30(7) 

Bohr-Weisskopf -0.0013(4) -0.061(18) -0.54(16) -8.2(2.5) -42(13) 

Quadrupole -0.0007(1) -0.018 -0.42 6.9(0.1) 7 

Recoil -0.0120 -0.015 -0.02 0.0 

Total 143.2960(5) 375.763(24) 656.36(18) 1456.3(2.6) 4046(15) 



The basic angular integrals Ai needed for the evaluation of the Zeeman vertex contribution are given by 

Ax{n a ,K h ) = J dpxi al/2 (p)voX K „i/2(p) , (A6) 

A 2 (K a ,K b ) = J dpx\ al/2 {f)i [<? x pioX Kl i/2(p) , (A7) 

A 3 (n a ,K b ) = J dpx[ al/2 (fi)i [a x Vn]oX K „ y 2 (p) , (A8) 

A 4 (K a ,K b ) = / dpx[ al/2 (f) i [p x Vn]ox Kb i/ 2 (f) ■ (A9) 



These integrals are evaluated by using the standard Racah algebra. The results relevant for this work are 

Ai(-l,-l) = l, ^(1,1) = -1/3, Ai(-1,2) = 0, Ai(l,-2) = -2v^/3, 

A 2 (-l, 1) = 2/3, A 2 (l, -1) = -2/3, A 2 (-l, -2) = V2/3, A 2 (l, 2) = -y/2/Z, 

A 3 (-l, 1) = 4/3, A 3 (l, -1) = 0, A 3 (-l, -2) = 2^2/3, A 3 (l, 2) = -A 

A 4 (-1,-1) = 0, Ai(l,l) = -2/3, A 4 (-l,2)=0, Ai(l,-2) = ->/2/3. 
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